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The Hubble relation between distance and redshift is a purely cosmographic 
relation that depends only on the symmetries of a FLRW spacetime, but does 
not intrinsically make any dynamical assumptions. This suggests that it should 
be possible to estimate the parameters defining the Hubble relation without 
making any dynamical assumptions. To test this idea, we perform a number of 
inter-related cosmographic fits to the legacyOS and gold06 supernova datasets, 
paying careful attention to the systematic uncertainties. Based on this super- 
nova data, the "preponderance of evidence" certainly suggests an accelerating 
universe. However we would argue that (unless one uses additional dynamical 
and observational information, and makes additional theoretical assumptions) 
this conclusion is not currently supported "beyond reasonable doubt" . As part 
of the analysis we develop two particularly transparent graphical representa- 
tions of the redshift-distance relation — representations in which acceleration 
versus deceleration reduces to the question of whether the relevant graph slopes 
up or down. 

Keywords: Cosmography, Hubble parameter, deceleration parameter, jerk. 



1. Introduction 

When analyzing the case for "dark energy" , it is critically important to 
realize that the standard luminosity distance versus redshift relation, 
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arc purely cosmographic results applicable to any FLRW universe, regard- 
less of the assumed dynamics. Following the spirit of Hubble's original pro- 
posal/ one could in principle fit such a relation directly to the supernova 
data,®"-^^ thereby estimating cosmological parameters (such as Hq, qq, and 
the jerk jo) without making any dynamical assumptions — but we shall 
see that there are ways of pre-processing the Hubble relation to make the 
result (and potential problems) stand out in greater clarity. 

For instance, it is sometimes obscrvationally more convenient to count 
photons rather than energy, and consider the "photon flux distance" ^''^^^ 



(1 + 0)V2' 



for which, defining 



dH = -^] and ^0 = 1 + = 1 + (4) 



one derives 

M^) = Ji - f + + ^ - + + 0(.')]. (5) 

Furthermore, using the "distance modulus" in terms of which the supernova 
data is actually reported^"^^ 

liD = 5 logio[dL/(10 pc)] = 5 logio[Q!L/(l Mpc)] + 25, (6) 

one has the simple relation 

ln[dF/iz Mpc)] = i^[/iD - 25] - Inz - i ln(l + z), (7) 
leading to a particularly useful form of the Hubble relation: 
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(8) 

Note that the question of whether or not the universe is accelerating or 
decelerating now reduces to the simple question of whether or not the curve 
obtained by plotting lB.[dF/{z Mpc)] versus z slopes up or down. 

For the data plots presented in this article we have used data from the 
supernova legacy survey (Iegacy05)^'^ and the Riess et. al. "gold" dataset 
of 2006 (gold06).^^ (The gold06 dataset is a larger dataset that contains 
most but not all of the Iegacy05 supernovae.) Note that figures 1 and 2 
are not plots of "statistical residuals" obtained after curve fitting — rather 
they can be interpreted as plots of "theoretical residuals" , obtained by first 
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Logarithmic Photon flux distance versus z-redshift using legacyOS 
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Fig. 1. The normalized logarithm of the photon flux distance, ln{dp/[z Mpc]), as a 
function of the 2-redshift using the Iegacy05 dataset.*'^ As is traditional in the field, the 
plotted error bars do not include estimates of the systematic errors. 



Logarithmic Photon flux distance versus z-redshift using gold06 
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Fig. 2. The normalized logarithm of the photon flux distance, ln{dp/[z Mpc]), as a 
function of the 2-redshift using the gold06 dataset.^"'^^ As is traditional in the field, the 
plotted error bars do not include estimates of the systematic errors. 
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splitting off the linear part of the Hubble law (which is now encoded in the 
intercept with the vertical axis), and secondly choosing the quantity to be 
plotted so as to make the slope of the curve at redshift zero particularly 
easy to interpret in terms of the deceleration parameter. 

The plots presented in figures 1 and 2 are considerably more ambiguous 
than we had initially expected. In generating these plots and performing the 
statistical analysis to be described below (considerably more detail can be 
found at"'^^"'^^) we had initially hoped to verify the robustness of the Hubble 
relation, and to possibly obtain improved estimates of cosmological param- 
eters such as the deceleration and jerk parameters, thereby complementing 
other recent cosmographic and cosmokinetic analyses such as,^^"^'^ as well 
as other analyses that take a sometimes skeptical view of the totality of the 
observational data.^^"^^ 

In view of the rather disturbing visual impact of figures 1 and 2 we 
resolved to see if they could be improved by further transformations of 
the data. For instance, we looked at the possibility of transforming the 
redshift variable, we looked at the possibility of adopting a number of other 
distance surrogates, and we performed a detailed statistical analysis of the 
data paying careful attention to the question of estimating the systematic 
uncertainties.^'* While the "preponderance of evidence" certainly suggests 
an accelerating universe, we would argue that (unless one uses additional 
dynamical and observational information, and makes additional theoretical 
assumptions) this conclusion is not currently supported "beyond reasonable 
doubt" . The supernova data (considered in isolation) certainly suggests an 
accelerating universe, but it is not suflacient to allow us to reliably conclude 
that the universe is accelerating. 



2. New redshift variable: y = z/(l + z) 

Because much of the recent supernova data is being acquired at large red- 
shift (z > 1), there are a number of theoretical reasons why it might be 
more appropriate to adopt the modified redshift variable^^"^^ 



''Prom recent discussions with a broad cross section of the community, it appears that this 
result now seems to have become part of the standard "folklore" . Typically, statistically 
strong arguments for cosmic acceleration rely on working within a particular dynami- 
cal framework (such as ACDM), and on extra observational data (such as independent 
constraints on Oq coming from CMB observations). 
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In the past (of an expanding universe) 

0G(O,oo); ye (0,1); (10) 

while in the future 

^€(-1,0); ye (-00,0). (11) 

Thus the variable y is both easy to compute, and when extrapolating back 
to the Big Bang has a nice finite range (0, 1). Furthermore, Taylor series 
in terms of the y variable have improved convergence properties at high 
rcdshift.^''"^^ We will refer to this variable as the y-redshift. ^ In terms of 
the variable y: 

. , ^ . /i [-3 + go] , [12 - 5go + 3gg - (jo + no)] 
dhiv) = dHy<l- - — - — -y + - 

It is now useful to define a quantity 

dh dp 

which we shall refer to as the "deceleration distance". This quantity has 
the nice feature that 

ln[rfQ/(y Mpc)] = - 25] - Iny + ^ ln(l - y), (14) 




whence 
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Thus plotting ln[dQ/(y Mpc)] versus y results in a curve whose slope at 
rcdshift y = is directly proportional to the deceleration parameter: The 
question of whether or not the universe is accelerating or decelerating now 
reduces to the simple question of whether or not the curve obtained by 
plotting ln[dQ/{y Mpc)] versus y slopes up or down. 

Visually, the plots presented in figures 3 and 4 arc again considerably 
more ambiguous than we had initially expected. Note that up to this point 
we have not performed any statistical analyses, we have "merely" found a 
dramatic way of visually presenting the observational data. 



''Similar expansion variables have certainly been considered before. See, for example, 

Chevalier and Polarski,^" who effectively worked with the dimcnsionlcss quantity b = 
a(t) / ao , so that y = 1 — b. Similar ideas have also appeared in several related works. ^''^^^ 
Note that these authors have typically been iut (-rested in parameterizing the so-called 
M!-parameter, rather than specifically addressing the Hubble relation. 
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Logarithmic Deceleration distance versus y-redsiiift using legaoyOS 




y-redstiift 



Fig. 3. The normalized logarithm of the deceleration distance, ln{dQ/[y Mpc]), as a 
function of the j;-redshift using the Iegacy05 dataset.*'® As is traditional in the field, the 
plotted error bars do not include estimates of the systematic errors. 



Logarithmic Deceleration distance versus y-redshift using gold06 
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Fig. 4. The normalized logarithm of the deceleration distance, ln(dQ/[y Mpc]), as a 
function of the y-rcdshift using the gold06 dataset.^'^'^^ As is traditional in the field, the 
plotted error bars do not include estimates of the systematic errors. 
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3. Data fitting: Statistical tests 

In view of the somewhat ambiguous and possibly alarming nature of the 
plots presented in figures 1-4, we then performed a number of statistical 
tests and analyses to check the extent to which robust and reliable results 
could be obtained from the supernova data considered in isolation. (For 
extensive technical details see/'*'^^ and related theoretical discussion in-^^). 
In performing the statistical analyses reported below we compared and 
contrasted results using several notions of cosmological distance, and two 
different versions of redshift. The distance surrogates we used were:^^"^^ 



• The "luminosity distance" cIl- 

• The "photon flux distance": dp = rfi (1 + z)-^^^ . 

• The "photon count distance": dp = dL {I + z)^^ . 

• The "deceleration distance": dg = di, (1 + z)~^/'^. 

• The "angular diameter distance": dA = di, (1 + z)~'^. 

The z-based versions of the Hubble law we used were:^^'"^^ 
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Similarly, the y-hased versions of the Hubble law we used were:^^'^^ 
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Fits were carried out for all five distance surrogates, and for both definitions 
of redshift, using polynomial approximants to the Hubble relation up to 7th- 
order. ^ The F-test was then used to discard statistically meaningless terms, 
and it was seen that quadratic fits were the best that could meaningfully 
be adopted. The results are presented in tables 1 4. 

Even after we have extracted these numerical results there is still a con- 
siderable amount of interpretation that has to go into understanding their 



•^Note that because the uncertainty in the redshift is encoded in the uncertainty of the 
distance modulus, the uncertainty in logarithmic distance is just scaled by a factor of 
ln{10)/5. Therefore, if the uncertainty is gaussian in the distance modulus, it is also 
gaussian in logarithmic distance, which is crucial for least squares fitting. 

'^Notc that in a cosmographic framework, where one is most closely following the spirit 
of Hubble's original methodology,'^ one does not have a dynamical model to fit the data 
to, and the use of least-squares fits to a truncated Taylor series is the best one can 
possibly hope for. Ultimately, the truncated Taylor series method is not really a very 
radical approach, being firmly based in quite standard statistical techniques.^''^^^ 
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Deceleration and jerk parameters (Iegacy05 dataset, y-redshift). 



distance 


<la 


jo + ^0 


dh 


-0.47 ±0.38 


-0.48 ±3.53 


dp 


-0.57 ±0.38 


+1.04 ±3.71 


dp 


-0.66 ±0.38 


+2.61 ±3.88 


dq 


-0.76 ±0.38 


+4.22 + 4.04 


dA 


-0.85 ±0.38 


+5.88 + 4.20 



With l-fT statistical uncertainties. 



Deceleration and jerk parameters (Iegacy05 dataset, z-redshift). 



distance 


Qo 


jo + ^0 


dp 


-0.48 + 0.17 


+0.43 + 0.60 


dp 


-0.56 + 0.17 


+1.16 + 0.65 


dp 


-0.62 + 0.17 


+1.92 + 0.69 


dQ 


-0.69 + 0.17 


+2.69 + 0.74 


dA 


-0.75 + 0.17 


+3.49 + 0.79 



With l-fT statistical uncertainties. 



Deceleration and jerk parameters (gold06 dataset, y-redshift). 



distance 


Qo 


jo + ^0 


dp 


-0.62 + 0.29 


+1.66 + 2.60 


dp 


-0.78 + 0.29 


+3.95 + 2.80 


dp 


-0.94 + 0.29 


+6.35 + 3.00 


dq 


-1.09 + 0.29 


+8.87 + 3.20 


dA 


-1.25 + 0.29 


+11.5 + 3.41 



With l-fT statistical uncertainties. 



Deceleration and jerk parameters (gold06 dataset, ^-redshift). 



distance 


Qo 


jo + ^0 


dp 


-0.37 + 0.11 


+0.26 + 0.20 


dp 


-0.48 + 0.11 


+1.10 + 0.24 


dp 


-0.58 + 0.11 


+1.98 + 0.29 


dq 


-0.68 + 0.11 


+2.92 + 0.37 


dA 


-0.79 + 0.11 


+3.90 + 0.39 



With 1-cr statistical uncertainties. 
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physical implications. In particular note that the differences between 
the various models, (Which distance do we use? Which version of redshift 
do we use? Which dataset do we use?), often dwarf the statistical uncer- 
tainties within any particular model. If better quality (smaller scatter) data 
were to become available, then one could hope that the cubic term would 
survive the i^-test. This would have follow-on effects in terms of making 
the differences between the various estimates of the deceleration parameter 
smaller, ^"''^^ which would give us greater confidence in the reliability and 
robustness of the conclusions. 

The fact that there are such large differences between the cosmological 
parameters deduced from the different models based on physically plausible 
distance indicators should give one pause for concern. These differences 
do not arise from any statistical flaw in the analysis, nor do they in any 
sense represent any "systematic" error, rather they are an intrinsic side- 
effect of what it means to do a least-squares fit — to a finite-polynomial 
approximate Taylor series — in a situation where it is physically unclear as 
to which if any particular measure of "distance" is physically preferable, and 
which particular notion of "distance" should be fed into the least-squares 
algorithm. (This "feature" — some may call it a "limitation" — of the least- 
squares algorithm in the absence of a clear physically motivated dynamical 
model is an often overlooked confounding factor in data analysis. ^^~^®) 



4. Systematic uncertainties 

In addition to the purely statistical uncertainties discussed above, one needs 
to make an estimate of the systematic uncertainties, and following NIST 
guidelines, *° combine the statistical and systematic uncertainties in quadra- 
ture 



Estimating systematic uncertainties is notoriously difficult. A careful de- 
scription of our own preferred way of estimating systematic uncertainties 
is fully discussed in,^**'^^ wherein we consider both modelling and histori- 
cal uncertainties. It should be emphasized that the (to our minds) overly 
optimistic estimates of systematic uncertainties commonly found in the lit- 
erature do not greatly change our conclusions below. 




(26) 
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5. Expanded uncertainties 

After due allowance is made for estimating the systematic imcertaintics, 
the NIST guidelines''^ recommend defining an "expanded uncertainty" by 

Uk = k CTcombincd- (27) 

Here the factor k is chosen for scientific (or legal) reasons to be such that 
one is "certain" that the true result lies within the stated range. The tra- 
dition within the social and medical sciences is to accept k = 2 (that is, 
two-sigma, corresponding approximately to 95% confidence intervals) as 
being sufficient to draw valid conclusions. Particle physics has traditionally 
adopted A; = 3 as the minimum standard for claiming "evidence for" a given 
hypothesis. (This is the origin of the aphorism: "If it's not threc-sigma, it's 
not physics".) Over the last 20 years or so, particle physics has moved to 
the more conservative consensus that fc = 5 is the minimum standard for 
claiming "discovery" of "new physics" . Our best estimates for the combined 
and expanded uncertainties are presented in tables 5-6. 



Deceleration parameter summary: Combined and expanded uncertainties. 



datascl 


rtHlsliift 




111) ± f^i 


± Ur, 


legacyOS 


y 


-0.66 ±0.42 


-0.66 ± 1.26 


-0.66 ±2.10 


Iegacy05 


z 


-0.62 ±0.23 


-0.62 ±0.70 


-0.62 ± 1.15 


goldOe 


y 


-0.94 ±0.39 


-0.94 ±1.16 


-0.94 ±1.95 


goldOe 


z 


-0.58 ±0.23 


-0.58 ±0.68 


-0.58 ± 1.15 



.Terk parameter suuuuary: Conibintxl and expanded uncertainties. 



dataset 


redshift 


{jo + f^o) ± O"combincd 


ijo + Oo) ± Us 


(jo + f^o) ± U5 


legacyOS 


y 


+2.65 ±4.63 


+2.65+13.9 


+2.65 ±23.2 


Iegacy05 


z 


+ 1.94 ±1.72 


+1.94 + 5.17 


+ 1.94 + 8.60 


goidoe 


y 


+6.47 + 4.75 


+6.47+14.2 


+6.47+23.8 


goidoe 


z 


+2.03+1.75 


+2.03 + 5.26 


+2.03 + 8.75 



6. Conclusions 

What can we conclude from this? While the "preponderance of evidence" 
is certainly that the universe is currently accelerating, qq < 0, this is not 
yet a "gold plated" result, at least not without bringing in other physical 
assumptions and observations; such as a specific dynamical model [e.g., 
ACDM] and/or invoking knowledge of flm or the CMB data, all of which 
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arc subject to their own additional theoretical assumptions. It is certainly 
more likely that the expansion of the universe is accelerating, than that 
the expansion of the universe is decelerating — but this is not the same as 
having definite evidence in favour of acceleration. 

We wish to emphasize the point that, regardless of one 's views on 
how to combine formal estimates of uncertainty, the very fact that 

different distance scales yield data-fits with such widely discrepant 
estimates for the cosmological parameters strongly suggests the need 
for extreme caution in interpreting the supernova data. 

There are a number of other more sophisticated statistical methods that 
might be applied to the data to possibly improve the statistical situation. 
For instance, ridge regression, robust regression, and the use of orthogonal 
polynomials and "loess curves" could all be adopted and adapted to focus 
more carefully on the region near rcdshift zcro.^^^'^^ However one should 
always keep in mind the difference between accuracy and precision. '^^ More 
sophisticated statistical analyses may permit one to improve the precision 
of the analysis, but unless one can further constrain the systematic un- 
certainties such precise results will be no more accurate than the current 
situation. 

However, we are certainly not claiming that all is grim on the cosmolog- 
ical front and do not wish our views to be misinterpreted in this regard 
— there are clearly parts of cosmology where there is plenty of high-quality 
data, and more coming in, constraining and helping refine our models. But 
regarding some specific cosmological questions the catch cry should still 
be "Precision cosmology? Not just yet" .'^^ In closing, we strongly encour- 
age readers to carefully contemplate figures 1-4 as an inoculation against 
over-interpretation of the supernova data. 

Ultimately, it is the fact that figures 1-4 do not exhibit any over- 
whelmingly obvious trend that makes it so difficult to make a robust 
and reliable estimate of the sign of the deceleration parameter. 

Finally we remind the reader that it is the putative acceleration of 
the expansion of the universe, no matter how derived, that then (via the 
Freidmann equations), is taken to imply the existence of "dark energy". 
In the absence of truly compelling model-independent evidence for cosmic 
acceleration one has to be at least a little cautious regarding the existence 
of "dark energy" . 
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